Two-resonator circuit QED: A superconducting quantum switch 
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We introduce a systematic formalism for two-resonator circuit QED, where two on-chip microwave 
resonators are simultaneously coupled to one superconducting qubit. Within this framework, we 
demonstrate that the qubit can function as a quantum switch between the two resonators, which are 
assumed to be originally independent. In this three-circuit network, the qubit mediates a geomet- 
ric second-order circuit interaction between the otherwise decoupled resonators. In the dispersive 
regime, it also gives rise to a dynamic second-order perturbative interaction. The geometric and 
dynamic coupling strengths can be tuned to be equal, thus permitting to switch on and off the 
interaction between the two resonators via a qubit population inversion or a shifting of the qubit 
operation point. We also show that our quantum switch represents a flexible architecture for the 
manipulation and generation of nonclassical microwave field states as well as the creation of con- 
trolled multipartite entanglement in circuit QED. In addition, we clarify the role played by the 
geometric interaction, which constitutes a fundamental property characteristic of superconducting 
quantum circuits without counterpart in quantum-optical systems. We develop a detailed theory of 
the geometric second-order coupling by means of circuit transformations for superconducting charge 
and flux qubits. Furthermore, we show the robustness of the quantum switch operation with respect 
to decoherence mechanisms. Finally, we propose a realistic design for a two-resonator circuit QED 
setup based on a flux qubit and estimate all the related parameters. In this manner, we show that 
this setup can be used to implement a superconducting quantum switch with available technology. 

PACS numbers: 03.67.Lx,42.50.Pq,84.30.Bv,32.60.+i 



I. INTRODUCTION 

In the past few years, we have witnessed a tremen- 
dous experimental progress in the flourishing realm of 
circuit QED.— i2i2il There, different types of supercon- 
ducting qubits have been strongly coupled to on-chip mi- 
crowave resonators, which act as quantized cavities. Re- 
cently, a quantum state has been stored and coherently 
transferred between two superconducting phase qubits 
via a microwave resonator— and two transition qubits 
have been coupled utilizing an on-chip cavity as a quan- 
tum bus£ Furthermore, microwave single photons have 
been generated by spontaneous emission 7 - and Fock states 
created in a system based on a phase qubit £ In addi- 
tion, lasing effects have been demonstrated exploiting a 
single Cooper-pair box^ the nonlinear response of the 
JC model observed , 10 ' 11 the two-photon driven Jaynes- 
Cummings (JC) dynamics used as a means to probe the 
symmetry properties of a flux qubit ^ and resonators 
tuned with high fidelity. 13 ' 14 These formidable advances 
show how circuit QED systems are rapidly reaching a 
level of complexity comparable to that of the already 
well-established quantum optical cavity QED. 15 ' 16 : 17 
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Amongst the aims common to these experiments is 
the possibility to perform quantum information process- 
ing,— in particular following the lines of recent propos- 
als, e.g., see Refs. [l^ and [2(J The latter considers a 
two-dimensional array of on-chip resonators coupled to 
qubits. In this or any other multi-cavity setup^ it is 
highly desirable to switch on and off an interaction be- 
tween two resonators or to compensate their spurious 
crosstalk. Moreover, investigating the basic properties of 
two- resonator circuit QED, where two resonators are cou- 
pled to one qubit, certainly represents a subject of fun- 
damental relevance. In fact, when operating such a sys- 
tem in a regime dominated by second-order (dispersive) 
interactions, as in the scope of this article, the require- 
ments on the qubit coherence properties are considerably 
relaxed. 22 ' 23 In this manner, two-resonator architectures 
constitute an appealing playground for testing quantum 
mechanics on a chip. We also notice that second-order 
interaction s) 12 ' 24 are becoming more and more prominent 
in circuit QED experiments owing to the possibility of 
very large first-order coupling strengths . 3 ' ' ' 9 ' ' i 

In this article, we theoretically study a three- 
circuit network where a superconducting charge or flux 
qubi t 26 ' 27 ' 28 ' 29 interacts with two on-chip microwave cav- 
ities, a two-resonator circuit QED setup. In the absence 
of the qubit, the resonators are assumed to have negligi- 
ble or small geometric first-order (direct) crosstalk. This 



2 



scenario is similar to that of quantum optics, where an 
atom can interact with two orthogonal cavity modest 
However, there are some crucial differences. The nature 
of the three-circuit system considered here requires to 
account for a geometric second-order circuit interaction 
between the two resonators. This gives rise to coupling 
terms in the interaction Hamiltonian, which are formally 
equivalent to those describing a beam splitter. This in- 
teraction is mediated by the circuit part of the qubit and 
does not depend on the qubit state. It is noteworthy 
to mention that this coupling does not exist in the two- 
mode Jaynes-Cummings (JC) model studied in quantum 
optics, where atoms do not sustain any geometric inter- 
action. This means that introducing a second resonator 
causes a departure from the neat analogy between cavity 
and one- resonator circuit QEDi 22 i 31 i 32 i 33 i 34 In the disper- 
sive regime, where the transition frequency of the qubit 
is largely detuned from that of the cavities, also other 
beam-splitter-type interaction terms between the two res- 
onators appear. Their existence is known in quantum 
optics^ and results from a dynamic second-order per- 
turbative interaction, which depends on the state of the 
qubit. The sign of this interaction can be changed by 
an inversion of the qubit population or by shifting the 
qubit operation point. The latter mechanism can also be 
used to change the interaction strength. Notably, for a 
suitable set of parameters, the geometric and dynamic 
second-order coupling coefficients can be made exactly 
equal by choosing a proper qubit-resonator detuning. In 
this case, the interaction between the two cavities can 
be switched on and off, thereby enabling the implemen- 
tation of a discrete quantum switch as well as a tunable 
coupler. 

In circuit QED, several other scenarios have been en- 
visioned where a qubit interacts with different bosonic 
modes, e.g., those of an adjacent nanomechanical res- 
onator or similar. It has been proposed to implement 
quantum transducers^ as well as Jahn- Teller models and 
Kerr nonlinearitiesj^l to generate nontrivial nonclassical 
states of the microwave radiation ) 21 i 38 to create entangle- 
ment via Landau-Zener sweeps^ and to carry out high 
fidelity measurements of microwave quantum fieldsi 38 ! 40 
Moreover, multi-resonator setups might serve to probe 
quantum walks 41 - and to study the scattering process 
of single microwave photons 42 All of these proposals, 
however, do not develop a systematic theory of a real- 
istic architecture based on two on-chip microwave res- 
onators and do take into account the fundamental geo- 
metric second-order coupling between them. Also, our 
quantum switch is inherently different from the quantum 
switches investigated in atomic systems 4 s - First, we con- 
sider a qubit simultaneously coupled to two resonators, 
which are not positioned one after the other in a cas- 
cade configuration as in Ref. 0. Second, our switch 
behaves as a tunable quantum coupler between the two 
resonators. Last, atomic systems naturally lack a ge- 
ometric second-order coupling. Furthermore, it is im- 
portant to stress that the dynamic interaction studied 



here cannot be cast within the framework of the quantum 
reactance theory (capacitance or inductance, depending 
on the specific implementation) 44,45,46,47,48,49,50,51 The 

main hypothesis for a quantum reactance to be defined 
is a resonator characterized by a transition frequency ex- 
tremely different from that of the qubit. Typically, the 
resonator frequency is considered to be very low (practi- 
cally zero) compared to the qubit one. Such a scenario 
is undesirable for the purposes of this work, where a 
truly quantized high-frequency cavity initialized in the 
vacuum state has to be used. Also, to our knowledge, 
the quantum reactance works mentioned above do not 
directly exploit a geometric coupling between two res- 
onators to compensate a dynamic one. Nevertheless, we 
believe that a circuit theory approach^Z.i^iSi^i^iS&S to 
two- resonator circuit QED, which we pursue throughout 
this manuscript, allows for a deep comprehension of the 
matter discussed here. Finally, we point out that the ge- 
ometric first-order coupling between two resonators can 
be reduced or erased by simple engineering, whereas the 
second-order coupling due to the presence of a qubit cir- 
cuit is a fundamental issue. As we show later, whenever 
the coupling between qubit and resonators is wanted to 
be large, an appreciable geometric second-order coupling 
inevitably appears, especially for resonators perfectly iso- 
lated in first order. In summary, our quantum switch is 
based on a combination of geometric and dynamic inter- 
actions competing against each other and constitutes a 
promising candidate to perform nontrivial quantum op- 
erations between different resonators. 

The paper is organized as follows. In Sec. [Til we 
develop a systematic formalism for two-resonator cir- 
cuit QED employing second-order circuit theory. In 
Sec. IIII1 we focus on the dispersive regime of two- 
resonator circuit QED and derive the quantum switch 
Hamiltonian. In Sec. HVi we discuss the main limitations 
to the quantum switch operation due to decoherence pro- 
cesses of qubit and cavities. In Sec. [V] we propose a 
realistic implementation of a two-resonator circuit QED 
architecture, which is suitable for the realization of a su- 
perconducting quantum switch. Finally, in Sec. IVI1 we 
summarize our main results, draw our conclusions, and 
give a brief outlook. 



II. TWO-RESONATOR CIRCUIT QED 

In this section, we take the perspective of classical cir- 
cuit theory^ and extend it to the quantum regime to 
derive the Hamiltonian of a quantized three-circuit net- 
work. In general, the latter is composed of two on-chip 
microwave resonators and a superconducting qubit. Our 
approach is similar to that of Refs. [271152.53 54,551. andl5rl 
In addition, we account for second-order circuit elements 
linking different parts of the network, which is consid- 
ered to be closed and nondissipativc. Here, closed means 
that we assume no energy flow between the network un- 
der analysis and other possible adjacent networks. These 
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FIG. 1: (Color online). Sketch of the system under analysis. All constants are defined in the main body of the paper. Only 
inductive couplings are considered, (a) Schematic representation of the first-order coupling Hamiltonian of our three-node 
network. Two cavities (resonators) A and B interact with a generic superconducting qubit Q. A and B can have a weak 
geometric first-order coupling M AB = m [broken blue (dark grey) arrow], as in the Hamiltonian H AB of Eq. (T2J. The two 
solid green (light grey) arrows represent a two-mode Jaynes-Cummings dynamics with coupling coefficients g A oc M A q and 
g B oc Mq B , respectively, (b) Visualization of the effective second-order coupling Hamiltonian H eS of Eq. (|14[) . The solid blue 
(dark grey) arrows show the second-order geometric coupling channel mediated by a virtual excitation of the circuit associated 
with Q, as in the Hamiltonian H AB of Eq. Q. This channel is characterized by a constant g AB oc M A qMq B /Mqq (the small 
contribution from m is neglected) and is qubit-state independent. The solid green (light grey) arrows show the second-order 
dynamic channel mediated by a virtual excitation of the qubit Q. This channel is characterized by a constant g A g B /A and is 
qubit-state dependent. Three generic sketches of a possible setup, (c) A flux qubit (Q) sits at the current antinode of, e.g., 
the first mode of two A/2 resonators (solid black lines, only inner conductor shown). The open circles at the "IN" and "OUT" 
ports denote the position of the coupling capacitors to be used in real implementations [e.g., cf. Sec. [V] and Fig. E^a)]. (d) A 
charge qubit (Q) sits at the voltage antinode of, e.g., the first mode of two A/2 resonators, (e) A charge or flux qubit sits at the 
voltage (e.g., second mode, A resonators) or current (e.g., first mode, A/2 resonators) antinode of two orthogonal resonators.— 



could be additional circuitry used to access the three- 
circuit network from outside and where excitations could 
possibly decay. Nondissipative means that we consider 
capacitive and inductive circuit elements only, more in 
general, reactive elements. We neglect resistors, which 
could represent dissipation processes of qubit and res- 
onators. In summary, the network of our model is alto- 
gether a conservative system. The detailed role of deco- 
herence mechanisms is studied later in Sec. IIVI 

The first step of our derivation is to demonstrate a 
geometric second-order coupling between the circuit ele- 
ments of a simple three-node network. This means that 
we assume the various circuit elements to be concentrated 
in three confined regions of space (nodes). Any topolog- 
ically complex three-circuit network can be reduced to 
such a three-node network, where each node is fully char- 
acterized by its capacitance matrix C and/or inductance 
matrix M. The topology of the different circuits (e.g., 
two microstrip or coplanar waveguide resonators coupled 



to a superconducting qubit) is thus absorbed in the defi- 
nition of C and M, simplifying the analysis significantly. 
The system Hamiltonian can then be straightforwardly 
obtained. In fact, the classical energy of a conservative 
network can be expressed as£= (V T CV + I T MI)/2, 
where the vectors V and I represent the voltages and 
currents on the various capacitors and inductors^ The 
usual quantization of voltages and currents^ and the ad- 
dition of the qubit Hamiltonian allows us to obtain the 
fully quantized Hamiltonian of the three-node network 
(cf. Subsec. Ill A[) . Special attention is then reserved to 
compute contributions to the matrices C and M up to 
second order. These arc consequently redefined as C' 2 - 1 
and M^ 2 ', respectively (cf. Subsec. fll B|) . Corrections of 
third or higher order to the capacitance and inductance 
matrices are discussed in Appendix f3 where we show 
that they are not relevant for this work. 

We finally consider two examples of possible implemen- 
tations of two-resonator circuit QED (cf. Subsec. fll C[) . 
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These examples account for two superconducting res- 
onators coupled to a charge quantum circuit (e.g., a 
Cooper-pair box or a transmon) or a flux quantum circuit 
(e.g., a superconducting one- or three- Josephson-junction 
loop). Before moving to a two- level approximation, the 
Hamiltonians of these devices can be used to deduce the 
geometric second-order circuit interaction between the 
two resonators. This result is better understood consid- 
ering the lumped-clement equivalent circuits of the en- 
tire systems. In this way, also the conceptual step from 
a three-circuit to a three-node network is clarified and 
the role played by the topology of the different circuits 
becomes more evident. We show that special care must 
be taken when quantizing the interaction Hamiltonian 
between charge or flux quantum circuits and microwave 
fields by the simple promotion of an AC classical field 
to a quantum one. Interestingly, comparing the stan- 
dard Hamiltonian of charge and flux quantum circuits 
coupled to quantized fields with ab initio models based 
on lumped-element equivalent circuits, we prove that the 
latter are better suited to describe circuit QED systems. 

A. The Hamiltonian of a generic three-node 
network 

The system to be studied is sketched in Figs. [TJa) and 
[ljb) , where the microwave resonators are represented by 
symbolic mirrors. A more realistic setup is discussed in 
Sec. |V] and is drawn in Fig. [6ja). A and B represent 
the two cavities and Q a superconducting qubit, making 
altogether a three-node network. The coupling channels 
between the three nodes are assumed to be capacitive 
and/or inductive. We also hypothesize the first-order in- 
teraction between A and B to be weak and that between 
A or B and Q to be strong by design. In other words, 
the first-order capacitance and inductance matrices are 
C = C u and M = M kl , with k, I G {A, B, Q}, where 
C kl = G lk and M kl = M lk because of symmetry rea- 
sons. In addition, we assume C AB = c -c G kl:jtAB and 
M AB = m <C Mki^AB- The elements c and m represent 
a first-order crosstalk between A and B, which can be ei- 
ther spurious or engineered and, here, is considered to be 
small. In Sec.[V] we delve into a more detailed analysis of 
the geometric first-order coupling between two microstrip 
resonators. Restricting the cavities to a single relevant 
mode, the total Hamiltonian of the system is given by 

h t = \v T cw v + X -P mW /+ \g (E c , e 3 ) b x , (1) 

where and M^"^ are the renormalized capacitance 
and inductance matrices up to the n-th order, with 

C' 1 ) ee C and Mt 1 ) ee M. Also, V = [V A ,V B ,V Q ] T and 

I ee [7 A , J B , Iq] t . In general, G is a function of the charg- 
ing energy E c and/or coupling energy Ej of the Joseph- 
son tunnel junctions in the qubit. For instance, G = 
for a charge qubit and G cx ^JE Q E^ cxp(—fj,y / E c /E- s ) 



for a flux qubit (/i ee const). Furthermore, V A = 

v dc + v ao( &1 + a)i = v Bofo + b )-. Vq = v q°z, 
[a = *dc + - a), / B ee i BO j(6 f - b), and 

Iq ee iq(7 2 . In these expressions, a x and a z are the 
usual Pauli operators for a spin- 1/2 system in the di- 
abatic basis, which consists of the eigenstatcs |— ) and 
1+) of CW'dc^q^ (charge case) or M AQ i BC i Q a z (flux 
case). Additionally, a) , w , a, and b are bosonic creation 
and annihilation operators for the fields of cavities A and 
B, respectively, and j = \/— 1. The DC voltage v BC and 
current i DC account for the quasi-static polarization of 
the qubit and can be applied through any suitable bias 
circuit. For definiteness, we have chosen here cavity A 
to perform this function. This is the standard approach 
followed by the charge qubit circuit QED community- 1 
However, for flux qubits the current i DC is more easily 
applied via an external coil.— i 12 i 25 i 58 In the latter case, we 
impose i DC = and add to the Hamiltonian of Eq. (fT]) 
the term ($° c — $ /2)Jq, where is an externally ap- 
plied flux bias and $ ee h/2e = 2.07 x 1CT 15 Wb is the 
flux quantum. The results of our derivation are not af- 
fected by this particular choice. The vacuum (zero point) 
fluctuations of the voltage and current of each resonator 
are give n by v A0 = x /fuj A /2C AA , v B0 = v / fuj B /2C BB , 

*ao - V^aT^aa, and *bo = V Huj b/ 2M bb> respec- 
tively. Here, lj a and lo b arc the transition angular fre- 
quencies of the two cavities. Finally, Vq and Iq represent 
the voltage of the superconducting island (s) and the cur- 
rent through the loop of the qubit circuit. Depending on 
the specific qubit implementation, either Vq or Iq domi- 
nates, thus defining the charge and flux regimes. 



B. The capacitance and inductance matrices up to 
second order 

The matrices C^ n ' and M' n ) account for corrections 
up to the n-th order interaction process between the ele- 
ments of the network. In fact, in order to write the exact 
Hamiltonian of the circuit, all possible electromagnetic 
paths connecting its nodes must be considered. A con- 
sequence of this approach to circuit theory is that the 
direct coupling 

H l Al=V A c% + I A mi B (2) 

between resonators A and B [cf. Fig. UK a)], here assumed 
to be small, is not the only interaction mechanism to be 
considered. In fact, an indirect coupling mediated by the 
circuit associated with the qubit Q has also to be included 
in the Hamiltonian. The dominating term for the A-Q-B 
excitation pathway can be derived from its second-order 
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electromagnetic energy [cf. Fig. QJb)], which gives 



f>(2) _ rr(l) 



Caq 7; Cqb Vb 

C>nn 



IaM aqj ^M qb I b . 



(3) 



Note that the inverse path (B-Q-A) is already in- 
cluded in this equation. In our work, we assume < 
c £ ^aqGqb/^qq and < m < M AQ M QB /M QQ 
(cf. Sec. |y]) . When c,m ~ 0, the direct coupling be- 
tween A and B is negligible, i.e., the contribution of H AB 
can be omitted. On the other hand, when c > and/or 
m > 0, both first- and second-order circuit theory con- 
tributions are relevant. In this case, c and m can rep- 
resent a spurious or an engineered crosstalk. The latter 
can deliberately be exploited to increase the strength of 
the geometric second-order coupling. However, c and m 
should be small enough to leave the mode structure and 
quality factors of A and B unaffected. 

From the knowledge of H AB , the capacitance matrix 
up to second order is readily obtained 



C (2) 
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(4) 



The second-order corrections to the self-capacitances, 
i.e., the diagonal elements C kk are absorbed in their def- 
initions^ (cf. Subsec. Ill C[) . In analogy, the corrected 
inductance matrix M' 2 ' is found substituting C kl with 
M kl and c with m in matrix ^ yielding 



second-order coupling, whose magnitude can be compa- 
rable to that of the geometric one, is also present (cf. Sub- 
sec. IIII Aj) . We study in more detail the relationship be- 
tween m and M A qMq B /Mqq in Sec. El There, we show 
that for a realistic design engineered for a flux qubit, 
which is our experimental expertise ; 12 ' 58 the geometric 
second-order interaction dominates over the first-order 
one. 

Figures Dlc)-(e) show three generic sketches, where the 
coupling of two on-chip resonators to one superconduct- 
ing qubit is illustrated. In particular, the sketch drawn 
in Fig. [He) is suitable when a flux qubit is intended 
to be utilized. In this case, the qubit is positioned at 
the current antinode of the first mode£°- of two A/2 res- 
onators. Moreover, this design clearly allows for engi- 
neering a strong coupling between the qubit and each 
resonator, while reducing the geometric first-order cou- 
pling between resonators A and B. This is due to the fact 
that the two cavities are close to each other only in the 
restricted region where the qubit is located and then de- 
velop abruptly towards opposite directions. The sketch 
in Fig. [Tf d) , instead, is more suitable for charge qubit 
applications. The qubit can easily be fabricated near a 
voltage antinodepii 2 ^ Similar arguments as in the previ- 
ous case apply for the qubit-resonator couplings and the 
geometric first-order coupling between A and B. Finally, 
the sketch of Fig. [TJe) relies on an orthogonal-cavity de- 
sign, which can be used for both charge and flux qubits. 
The main properties of such a setup have already been 
presented in two of our previous works ; 20 ' 21 where orthog- 
onal cavities have been exploited for different purposes. 
In conclusion, we want to stress that based on the gen- 
eral sketches of Figs, [ljc)-(e), a large variety of specific 
experimental implementations can be envisioned. 



C. The role of circuit topology: Two examples 
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(5) 

Again, second-order corrections to the self-inductances 
are absorbed in the definition of M kk . The matrices 
C' 2 - 1 and M' 2 ) constitute the first main result of this 
work. They show that, if a large qubit-resonator cou- 
pling (i.e., a vacuum Rabi coupling oc C A q , Cq B for 
charge quantum circuits and oc M A q , Mq B for flux 
quantum circuits) is present, as in most circuit QED 
implementations ±M£&L&2J2AL11 I 22& a relevant geo- 
metric second-order coupling (oc C A qCq B /Cqq or oc 
M A qMq B /Mqq for charge and flux quantum circuits, 
respectively) has to be expected. This coupling becomes 
relevant in the dispersive regime ; 22 ' 24 where a dynamic 



All results of Subsecs. Ill Al and III Bl are general and do 
not rely a priori on the knowledge of the three-circuit 
network topology. Here, we explain with the aid of two 
easy examples how to obtain a reduced three-node net- 
work starting from a three-circuit one. The examples are 
based on the coupling of two superconducting coplanar 
waveguide or microstrip resonators to a single Cooper- 
pair bosii 2 ^ (or a transmo n 61 ' 62 ' 63 ) or to a superconduct- 
ing loop interrupted by one (or three) Joscphson-tunnel 
junctions ^32^64 

The first example is the case of a single Cooper-pair 
box (a charge quantum circuit), which is formally equiva- 
lent to the more appealing case of the transmon. A single 
Cooper-pair box*^ is made of a superconducting island 
connected to a large reservoir via two Josephson tunnel 
junctions with Josephson energy Ej and capacitance Cj. 
The box is capacitevely coupled to two resonators A and 
B by the gate capacitors C ga and C gb , respectively. In 
the charge basis, the Hamiltonian of a single Cooper-pair 
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FIG. 2: Equivalent circuit diagrams for two different implementations of two-resonator circuit QED based on either a charge 
qubit [(a)-(c)] or a flux qubit [(d)-(f)]. Cf. Subsec. Ill CI for details, (a) V IlL and V rh : Quantized voltage sources associated 
with resonators A and B in parallel to the self-capacitances C ra and C rb of the resonators. 7 ra and I lb : Quantized current 
sources associated with resonators A and B in series to the self-inductances L ra and L rh of the resonators. The number of 
excess Cooper pairs on the charge qubit island (big dot) is (n\h\n). C 3 : Capacitance of each of the two Josephson tunnel 
junctions connecting the island to ground. C ga and C gb : Coupling capacitances between the qubit and the two resonators. C ab : 
First-order cross-capacitance between A and B (typically small, dotted line). The dashed box marks a T-network composed 
of C ga , 2Cj, and C gb . (b) C cr = C ga C gb /C E : Second-order cross-capacitance. C sa = 2CjC ga /C E and C sb = 2CjC gb /C s : 
Resonator shift capacitances. C ab is neglected for simplicity, (c) The circuits of (b) rearranged as a single EE-network (dashed 
box). The latter is equivalent to the T-network of (a). The magnitudes of C ra and C rb are increased by the presence of the 
shift capacitances C aa and C ab . (d) Two resonators A and B inductively coupled via M qa and M qb to a flux qubit with total 
self-inductance L q = L q /2 + L q /2 and flux operator $. The first-order mutual inductance m between the two resonators is 
neglected to simplify the notation, (e) The disconnected circuit of (d) is transformed into a connected circuit.— Again, we can 
identify a Fi-network (dashed box), (f) Left side: T-network obtained from the FT network of (e). We identify the second-order 

M qb /L q . Right side: The connected 



mutual inductance M cr = M qa Af qb /L q and the shift inductances L S3 
circuit on the left side is transformed into a disconnected circuit.— 



M qa /L q and L sb 



box can be written as=2 

H c = 4E c J2(n-n g ) 2 \n)(n\ 

n 

-^X>x n+1 i+i n+i x n i)' (6) 

n 

where E c = e 2 /2C s is the box electrostatic energy (e is 
the electron charge), C E = C ga + 2C 7 + C gb is its total 
capacitance^ (n| n \n) represents the number of excess 
Cooper pairs on the island, and n g is the global dimen- 
sionless gate charge applied to it. The latter is the sum of 
a DC signal c (here, considered to be applied through 
cavity A) and a high-frequency excitation (5n g applied 
through cavities A and/or B, n g = n^ c + Sn g . In par- 
ticular, Sn g can represent the quantized electric fields 
(equivalent to the voltages) of the two cavities acting as 
quantum harmonic oscillators. Restricting ourselves to 



the two lowest charge states n = 0, 1, we can rewrite the 
Hamiltonian of Eq. © as 

H c = 2E C (1 - 2n g + 2n 2 g + Z z - 2n g t z ) - ^ x 

= 2E C (1 - 2n° c ) b z - ^ x 

- 4E c Sn g (1 - 2n^ c - Sn g + & z ) . (7) 

The second line of the above equation forms the standard 
charge qubit, which can be controlled by the quasi-static 
bias ?ig C = C ga v DC /2e. The third line contains four 
high-frequency interaction terms. Among those, two of 
them are particularly interesting^ These are AE c Sn g 
and — AE c Sn g a z . We now quantize the high-frequency 
excitations Sn g — * 6n g = C ga w AO (a^-|-a)/2e+C gb w B g(6^-|- 

b)/2e, using the fact that they are the quantized volt- 
ages of the two resonators. We subsequently perform a 
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rotating-wave approximation (RWA) and, finally, write 
the interaction Hamiltonian 

= hG AB (tf + a)(P + b) 

- hG A a z (a ] +&)- hG B a z (P + b) 

+ huj A cra + hij B b'b , (8) 

where all constant energy offsets, e.g., the Lamb shifts, 
have been neglected. Remarkably, in the first line of 
the above equation we identify a geometric resonator- 
resonator interaction term with second-order coupling 
coefficient G AB = v A0 v B0 C g& C gh /C^h. Furthermore, 
the two terms of the second line of this equation repre- 
sent the expected first-order qubit-resonator interactions 
with coupling coefficients G A = e(C„ a /C s )v A0 /S and 
G B = e(C gb /C s )w B0 /?i, respectively. In the third line, 
the two small renormalizations uj a = {C gll v A0 ) 2 jC^h, and 
u5 B = (C gh v B0 ) 2 /C-^h of the resonator angular frequen- 
cies are artifacts due to the simple model behind the 
Hamiltonian of Eq. ([6]). A more advanced model based 
on a realistic circuit topology yields similar renormal- 
ization terms, which, however, are governed by different 
topology-dependent constants. Among the possible ways 
to find the correct constants, we choose the circuit trans- 
formations of Figs. [2](a)-(c). This approach also allows 
us to better understand the geometric second-order in- 
teraction term. 

In Fig. [2a), the two cavities are represented as LC- 
resonators with total capacitances and inductances C ra , 
C rb , L ra , and L b , respectively The quantized voltages 
and currents of the two resonators are V I& = v A0 (a ji + ®)> 
Kb = v bo(^ + b), l ra = i A0 j(al - a), and f rb = 
i B0 j(6' i ' — b), respectively Also, C ab accounts for a 
first-order cross-capacitance between resonators A and 
B, which, for simplicity, is neglected in Eqs. (J7]) and ([8]). 
In addition, here we are only interested in the geomet- 
ric properties of the charge quantum circuit. The dy- 
namic properties of this circuit are studied following a 
more canonical approach within a two-level approxima- 
tion in Sec. IIIII The dynamic properties arc governed 
by the two Josephson tunnel junctions and by the num- 
ber of excess Cooper pairs on the island, (n\n\n). To 
simplify our derivations, we can then assume n = and 
consider only the capacitance Cj of the two Josephson 
tunnel junctions, but not their Josephson energy. 

We now derive in three steps the geometric part of the 
interaction Hamiltonian by means of circuit theory. The 
procedure is visualized in Figs. (Ha)-(c). The steps are: 

(i) - First, we assume that the circuit associated to the 
charge qubit is positioned at a voltage antinodo^ 2 - of both 
resonators. Consequently, we can replace the two current 
sources of Fig. [^a) with open circuits, 7 ra = / rb = 0. 
Thus, we can eliminate both L ra and L rb from the circuit 
diagram because they are in series to open circuits. 

(ii) - Second, we apply the superposition principle of 
circuit theory^i One at the time, wc replace each of 
the two voltage sources with short circuits, V Ta = or 



V rh = 0. This allows us to split up the circuit of Fig. [2(a) 
into the two subcircuits of Fig. [DJb), which are topolog- 
ically less complex. As a consequence, in the respective 
subcircuits, C rb or C ra can be substituted by short cir- 
cuits and all other capacitors opportunely rearranged. 
In this way, for the case of cavity A, we find the small 
shift capacitance C sa = 2CjC ga /C s , which gives the cor- 
rect angular frequency renormalization of the resonator, 
~corr = 2CjC ga v A0 /C s h. Remarkably, we also find the 
second-order cross-capacitance C cr = C ga C gb /C E , corre- 
sponding to the geometric second-order coupling between 
the resonators. This coincides with our result obtained in 
Eq. ([3]) of Subsec. Ill Bl and is consistent with the simple 
model of Eqs. ©-©. We notice that C cr deviates from 
the simple series of the two gate capacitances C ga and C gb 
because of the presence of Cj in C s . For the case of cavity 
B, C sh ee 2CjC gb /Cx and Z B °" = 2C,C gh v 2 B0 /C s h can 
be derived in an analogous manner. In Subsec. Ill Bl the 
two renormalization constants as well as Cj are absorbed 
in the definitions of C AA , C BB , and Cqq, respectively. 

(Hi) - Third, we notice that the cross-capacitance C cr , 
which is responsible for the geometric second-order inter- 
action between A and B, is subjected to both quantum 
voltages V la and V rh . Hence, we can finally draw the 
circuit diagram of Fig.[^c). Indeed, we could have iden- 
tified the T-network of Fig. [DJa) (indicated by a dashed 
box) and transformed it into the equivalent n-network of 
Fig. [^c) (also indicated by a dashed box) in one single 
step^i obtaining the same results. We prefer to explicitly 
show the steps of Fig. [Hb) for pedagogical reasons. 

The second example is based on a superconducting 
loop interrupted by one Josephson tunnel junction (a 
flux quantum circuit). Such a device is also known as 
radio-frequency (RF) superconducting quantum interfer- 
ence device (SQUID). We choose the RF SQUID here 
for pure pedagogical reasons. In fact, our treatment 
could be extended to the more common case of three 
j unctions The Hamiltonian of an RF SQUID can be 
expressed a a 26 i 27 ' 32 i 60 

where Q is the operator for the charge accumulated on 
the capacitor Cj associated with the Josephson tunnel 
junction. The flux operator $ is the conjugated variable 
of Q, i.e., [$, Q] = jh. In analogy to the dimension- 
less gate charge n g of the previous example, the flux 
bias <f> x = <f>° c + <5(f> x consists of a DC and an AC 
component. The sclf-inductancc of the superconducting 
loop is defined as L q . When the RF SQUID is cou- 
pled to two quantized resonators, we can quantize the 
high-frequency excitations performing the transforma- 
tions 5$ x -> 5$ K = Af qa i A0 j(ot - o) + M qb i BO j(6+ - S). 
Here, M qa and M qb are the mutual inductances be- 
tween the loop and each resonator. We can then assume 
$ = 0, perform a two-level approximation and a RWA, 
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and finally obtain the same interaction Hamiltonian as 
in Eq. ©. However, in this case the coefficients are re- 
defined as Zj a ee (M q J AQ ) 2 /L q h, Zj b ee (M qh i m ) 2 / L q h, 
and G AB = i A0 i B0 M qa M b / L q h. The term with cou- 
pling coefficient G AB constitutes the geometric second- 
order interaction between A and B. As it appears clear 
from the discussion below, once again the renormaliza- 
tion terms uj a and lu b do not catch the circuit topology 
properly. This issue can be clarified analyzing the cir- 
cuit diagram drawn in Fig. EJd), where all the geometric 
elements for this example are shown. The geometric first- 
order mutual inductance m between the two resonators is 
neglected to simplify the notation. Again, the Josephson 
tunnel junctions responsible for the dynamic behaviour 
are not included. 

We now study the geometric part of the interaction 
Hamiltonian between the RF SQUID and the two res- 
onators following a similar path as for the case of the 
single Cooper-pair box [cf. Figs.[2jd)-(f)]. The four main 
transformation steps are: 

(i) - First, we assume the circuit corresponding to the 
flux qubit to be positioned at a current antinode. Thus, 
in Fig.[2jc), we replace all voltage sources and capacitors 
of Fig. HJd) with short circuits. The self-inductance of 
the qubit loop is split up into two L q /2 inductances to 
facilitate the following transformation steps. 

(ii) - Second, a well-known theorem of circuit theory^ 
allows us to transform the three disconnected circuits of 
Fig. [UJd) into the connected circuit of Fig. [He). Here, 
the region indicated by the dashed box evidently forms 
a n-network. 

(Hi) - Third, a n-to-T- network transformation^ re- 
sults in the circuit on the left side of Fig. [2jf ) . 

(iv) - Fourth, applying the inverse theorem of that 
used in step (ii) finally allows us to draw the equiva- 
lent circuit on the right side of Fig. [210 • Here, M CI = 
M qa M qh /L q represents the second-order mutual induc- 
tance between resonators A and B, corresponding to the 
geometric second-order coupling between them. Remark- 
ably, this coincides with our result obtained in Eq. ([3]) 
of Subsec. Ill Bl and is consistent with the simple model 
of Eq. ([9]). However, in this model the small shift in- 
ductances L sa ee M\jL q and L sh = M 2 h /L q (here de- 
fined to be strictly positive) acquire the wrong sign. Our 
circuit approach reveals that the correct renormaliza- 
tion constants of the resonators angular frequency are 
lo a = — L s& i 2 A0 /h and u> B = — L sh i B0 /h. This result is 
also confirmed by our numerical simulations (cf. Sec. IVl 
and Table [J). In Subsec. Til Bl these renormalization con- 
stants are absorbed in the definitions of M. , and M RR . 



III. DERIVATION OF THE QUANTUM 
SWITCH HAMILTONIAN 

In this section, we analyze the Hamiltonian of a three- 
node quantum network as found in Subsec. Ill Bl In par- 
ticular, we focus on the relevant case of large qubit- 



resonator detuning, i.e., the dispersive regime of two- 
resonator circuit QED. Under this assumption, we are 
able to derive an effective Hamiltonian describing a quan- 
tum switch between two resonators. Wc compare the an- 
alytical results to those of extensive simulations (cf. Sub- 
sec. IIII A|) . We also propose a protocol for the quantum 
switch operation stressing two possible variants. One is 
based on a qubit population inversion and the other on an 
adiabatic shift pulse with the qubit in the energy ground- 
state (cf. Subsec. IIII Bp . Finally, we give a few exam- 
ples of advanced applications of the quantum switch and, 
more in general, of dispersive two-resonator circuit QED 
(cf. Subsec. HHO]). 



A. Balancing the geometric and dynamic coupling 

We now give the total Hamiltonian of the three-node 
quantum network of Figs. [2(a) and [Hb). In order to 
avoid unnecessarily cumbersome calculations, we restrict 
ourselves to purely inductive interactions up to geomet- 
ric second-order corrections. In this framework, the most 
suitable quantum circuit to be used is a flux qubit. Here- 
after, all specific parameters and corresponding simula- 
tions refer to this particular case. Nevertheless, the for- 
malism which we develop remains general and can be ex- 
tended to purely capacitive interactions (charge qubits) 
straightforwardly. 

The flux qubit is assumed to be positioned at a current 
antinode. As a consequence, the vacuum fluctuations i A0 
and i B0 have maximum values »™q x and i B Q X at the qubit 
position and we can impose v A0 = v B0 = 0. Also, in 
the standard operation of a flux qubit no DC voltages 
are applied, i.e., i> DC = 0, and the quasi-static flux bias 
is usually controlled by an external coil and not by the 
cavities (cf. Subsec. Ill A [ I. Again, we impose i DC = 
and add to the Hamiltonian of Eq. (Q]) the term ($° C — 
$ /2)jq. Under all these assumptions and substituting 
M<") of Eq. (U) by M"( 2 ) of matrix ©, we readily obtain 

H' = —hea z + -Mqffj + hLo A a)a + ku! B b^b 

+ frgA&z(^ +&) + hg B a z (tf + b) 

+ hg AB {a) + a){P + b). (10) 

Here, all global energy offsets have been neglected and we 
have included both first- and second-order circuit the- 
ory contributions. In this equation, & ee 2ig($^ — 
$ /2) is the qubit energy bias, (5q ee 8q(E c ,Ej) is 
the qubit gap ) 26 ' 67 oj a = ^/ \/M AA C AA and uj b ee 
1/ 'y/ M BB C BB are the angular frequencies of resonators 
A and B, respectively, g A = «q*ao-^aq/^ anc ^ 5b = 
*q*bo^bq/^ are t ne qubit-resonator coupling coeffi- 
cients, and, finally, the second-order coupling coefficient 

Sab = Wbo ( m + M aq m qb/ m qq) l h - In general, g A 
and g B can be different due to parameter spread during 
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FIG. 3: (Color online) Simulation of the Hamiltonian of Eq. (|10[) in the dispersive regime (cf. Subsec. MI Al for a detailed 
description of the system parameters) . (a) The differences between the first nine excited energy levels of the quantum switch 
Hamiltonian and the groundstate energy level, AE, as a function of the frustration f^ c = $x C /<I> . The two lowest lines 
[blue (dark grey) and green (light grey), respectively] are associated with resonators A and B. The dispersive action of the 
qubit, which modifies the shape of the resonator lines, is clearly noticeable in the vicinity of the qubit degeneracy point. In 
this region, the third energy difference (hyperbolic shape, red line) represents the modified transition frequency of the qubit. 
(b) Close-up of the area indicated by the black arrow in (a). Here, the two modified resonator lines [thick blue (dark grey) and 
thin green (light grey), respectively] cross each other, (c) Quantum switch coupling coefficient |2p|w| extrapolated from the 
energy spectrum of (a) plotted versus /jP c . (d) Close-up of the area indicated by the black arrow of (c). The switch setting 
condition \2g l £ \ = is reached at f°° ~ 0.4938. 



the sample fabrication. Later, we show that the archi- 
tecture proposed here is robust with respect to such im- 
perfections. We now rotate the system Hamiltonian of 
Eq. (|10p into the qubit energy eigenbasis {|g) , |e)} ob- 
taining 

H = h—^-cr z + hjjpjv a + huj-ghrb 

+ Hg A cos 9a z (a' + a ) + hg B cos 9a z (b' + b) 
— Hg A sm9a x (a^ + a) — hg B cos9a x (W + b) 
+ % AB (a t + a)(& t +S). (11) 



Here, fiq = J e 2 + 8q is the $^ c -dependent transition 

frequency of the qubit and 9 = arctan ((5g/e) is the usual 
mixing angle. Hereafter, we use the redefined Pauli oper- 
ators a x and er 2 , where a x = a + +&~ and & + and &~ are 
the lowering and raising operators between the qubit en- 



ergy groundstate |g) and excited state |e), respectively. 
Expressing H in an interaction picture with respect to 
the qubit and both resonators, a) — > a) exp (+joj A t), a — * 
aexp(~juj A t), 6+ ->■ &t exp (+ jw B t) , b -> b exp (-ju B t), 



a T — > <t t exp (=Fjf2qf ) , assuming ui A 
and performing a RWA yields 



H = hsh 



9a 



a^+g B py ]At +a+(g A a+g B \ 



%ab (a f 



b+atf 



2tt/, 



(12) 



Here, A = J7q — w is the qubit-resonator detuning. The 
first two terms of Eq. (fT2"|) represent a standard two-mode 
JC dynamics i 30 i 69 The last term, instead, constitutes a 
beam-splitter-type interaction specific to two-resonator 
circuit QED. This interaction is not present in the quan- 
tum optical versioni 30 i 69 The coupling coefficient <7 AB is 
typically much smaller than g A and g B (see below) . How- 
ever, in the dispersive regime (|A| 3> max{(; A , g B , <?ab})> 
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Sab becomes comparable to all other dispersive coupling 
strengths. To gain further insight into this matter, we 

can define two superoperators = a + (g A a + g p b) and 



g B b^). It can be shown that the Dyson 



9a c 

series for the evolution operator associated with the time- 
dependent Hamiltonian of Eq. §12\i can be rewritten 



in the exponential form U = exp I —jH cS t/h\, where 



/A + hg AB [a)b + abU. Thus 



H oS = h 



H eS = h ^ ° z { a + 2 



• (13) 



In this Hamiltonian, the first two terms represent dy- 
namic (AC-Zeeman) shifts (AC-Stark shifts in the case 
of charge qubits) of the transition angular frequency of 
resonators A and B, respectively. If g A = g B = g and we 
only use eigenstatcs of <x z , the first two terms of Eq. (fTS")) 
equally renormalize uj a and lu b , respectively. The Hamil- 
tonian of Eq. (fT3"|) can be further simplified through an 
additional unitary transformation described by U = 



exp(jH t/H), where H = h(g A sxn 9/A)a z ( 



at; 



■1/2)- 



h(g B sm 2 6/A)a z (tfb + 1/2). When g A = g B = g, this 
transformation yields the final effective Hamiltonian 



h cS = n 



g 2 sin 2 i 



Z* + 9ab) (& r b + ab^ , (14) 



which constitutes the second main result of this work. 
This Hamiltonian is the key ingredient for the implemen- 
tation of a quantum switch between the two resonators. 
In fact, it clearly represents a tunable interaction between 
A and B characterized by an effective coupling coefficient 



, g> , 9 2 sin 2 i 

</sw — i/AB a 



o) _ 9 2 sin 2 I 

SW = g AB H - 



(15) 



for |g) and |e), respectively. The switching of such an in- 
teraction triggers, or prevents, the exchange of quantum 
information between A and B. On the one hand, the first 
part of this interaction is a purely geometric coupling, 
which is constant and qubit-state independent. On the 
other hand, the second part is a dynamic coupling, which 
depends on the state of the qubit. The switch setting 
condition 



g 2 sin 2 I 



I Sab I 



(16) 



can easily be fulfilled varying A, changing sin (9, or in- 
ducing AC-Zeeman or -Stark shifts.— In the special case 
of a charge qubit, not treated here in detail, this task 
can also be accomplished modifying the qubit transition 
frequency f2q via a suitable quasi-static magnetic fieldi 
This allows one to keep the qubit at the degeneracy point. 
Here, we focus on the first option, i.e., finding a suitable 
qubit bias for which the detuning A fulfills the relation of 
Eq. (|16p . For a flux qubit, this can be realized polarizing 
the qubit by means of an external flux. 

To better understand the switch setting condition, 
we numerically diagonalize the entire system Hamilto- 
nian of Eq. (|10p . without performing any approxima- 
tion. The results are presented in Fig. [3j which shows 
the energy spectrum of the quantum switch Hamilto- 
nian and the effective coupling coefficient \2gi^ | for a flux 
qubit with i Q = 370 nA, 5 Q /2n = 4 GHz, '/ = 3.5 GHz, 
g/2n = 472 MHz, and g AB /2n = 2.2 MHz. The parame- 
ters for the flux qubit are chosen from our previous ex- 
perimental works , 12 ' 58 whereas the three coupling coeffi- 
cients are the result of detailed simulations (cf. Sec. |V|) . 
It is noteworthy to mention that large vacuum Rabi cou- 
plings g/2ir on the order of 500 MHz have already been 
achieved both for flux and charge qubits.— ^ We have cho- 
sen the qubit to be already detuned from both resonators 
by 0.5 GHz when biased at the flux degeneracy point. 
Moving sufficiently far from the degeneracy point enables 
us to increase the qubit-rcsonator detuning such that the 
system can be modeled by the Hamiltonian of Eq. (|14[) . 
Figure 03a) shows the differences between the first nine 
excited energy levels of the quantum switch Hamilto- 
nian and the groundstate energy level, AE i = E i — E Q 
with i — {1, ...,10}, as a function of the frustration 
/ x = $x C /$ . Here, E i is the energy level of the i-th 
excited state and E that of the groundstate. Due to 
the qubit-resonator detuning, the two lowest energy dif- 
ferences [blue (dark grey) and green (light grey) lines, 
respectively] correspond to the modified transition fre- 
quencies of the two resonators. Owing to the interaction 
with the qubit these lines are not flat. This effect be- 
comes particularly evident in the region close to the qubit 
degeneracy point, where dispcrsivity is reduced. In this 
region, the third energy difference (hyperbolic shape, red 
line) represents the modified transition frequency of the 
qubit. When moving away from the qubit degeneracy 
point, a crossing between the modified resonator lines is 
encountered, as clearly shown in Fig.[3^b) [see, thick blue 
(dark grey) and thin green (light grey) lines] . This cross- 
ing represents the switch setting condition of Eq. (| 16|) . 
Figurcs[3]Jc) and[3][d) show the absolute value of the flux- 
dependent coupling coefficient \2gl^ \ in the flux windows 
of Figs. 03a) and[3{b), respectively. The switch setting 
condition \2g l £ \ = is reached at f^ c ~ 0.4938. 

A comparison between the analytic expression of 
Eq. (fT5")) with the qubit in |g) [dashed green (light grey) 
lines] and a numerical estimate of the effective cou- 
pling coefficient 1 2 w | [solid blue (dark grey) lines] is 
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FIG. 4: (Color online) Comparison between the /° c -dependence of the analytical expression for the coupling coefficient |2(?lw| 
obtained from Eq. (|15|l and the one found by numerically diagonalizing the Hamiltonian of Eq. (|10[) . (a) We choose a center 
frequency f A — f B = / = 2.7 GHz for the two resonators. All the other parameters are the same as those used to obtain 
the results of Fig. [3] The analytical [dashed green (light grey) line] and the numerical [solid blue (dark grey) line] results 
are in excellent agreement. In the large detuning limit far away from the qubit degeneracy point, |2griw| saturates to the 
value |2<7abI — 2.6 MHz. Inset: Close-up of the region near the switch setting condition, (b) Here, we choose a center 
frequency f A = / B = / = 3.5 GHz for the two resonators. The analytical [dashed green (light grey) line] and numerical 
[solid blue (dark grey) line] results are in good agreement away from the qubit degeneracy point. Closer to it they diverge 
(cf. Subsec. IlII Al for more details). In the large detuning limit far away from the qubit degeneracy point, |2p|w| saturates to 
the value 2<? AB | ~ 4.4MHz. Inset: Close-up of the region near the switch setting condition. 



shown in Figure 2] To clarify similarities and differ- 
ences between analytical and numerical calculations, we 
choose two different sets of parameters. In Fig. |U[a), 
the center frequencies of the two resonators are set to be 
/a = /b = / = 2.7 GHz, whereas in Fig. |3Jb) we choose 
f A = / B = / = 3.5 GHz. All the other parameters are 
equal to those used to obtain the results of Fig. [31 In 
Fig. SJa), analytics and numerics agree over the entire 
frustration window. The inset shows that the switch set- 
ting condition obtained from the numerical simulation is 
only slightly shifted with respect to the analytical predic- 
tion. Also in Fig.QJb), the agreement between analytical 
and numerical estimates is good far away from the qubit 
degeneracy point. However, closer to it the qubit and 
the two resonators are not detuned enough to guarantee 
dispcrsivity. Therefore, analytics and numerics start to 
deviate, as expected. Again, the inset shows that the 
switch setting condition can be fulfilled. It is noteworthy 
to point out that both analytical and numerical estimates 
converge to the value |2g AB | in the limit of large detuning. 
We find |2.g AB |/27r ~ 2.6MHz and |2 5ab |/2tt ~ 4.4MHz 
from the simulations that produce Figs. Ufa) and0|b), 
respectively. 

Finally, we demonstrate that the quantum switch 
Hamiltonian is robust to parameter spread due to fab- 
rication inaccuracies. Typically, for a center frequency 
of 5 GHz the expected spread around this value is ap- 
proximately^ =f10MHz for two resonators fabricated on 
the same chip. Also, the coupling coefficients g A and g B 



can differ slightly. In this case, a generalized effective 
Hamiltonian for the quantum switch can be derived 2 ^ 



f^gon fc (g A sin6>) „ . 



A 



A 



9a9b ! 



1 



Sab 



x (atg e +j«AB* + aPe~ jS ^ , (17) 



where A A = £l c 



-A- 



A R EE Sir 



j b , and <S AB ee 



lu a — lo b . From Eq. (jTTJ) , we can deduce the general- 
ized coupling coefficient of the switch, g!w'' G ^ = g AB T 
g A g B sin 2 9 (1/2 A A + 1/2 A B ) for the qubit groundstate 
|g) or excited state |e), respectively. As a consequence, 
the generalized switch setting condition becomes 



9a9b sm 



1 



1 



= IffABl 



(18) 



This condition is displayed in Fig. [5] [dashed green (light 
grey) line] as a function of the external flux. Here, 
we assume two resonators with center frequencies f A = 
3.5 GHz and f B = 3.5 GHz + 35 MHz. This corre- 
sponds to a relatively large center frequency spread^ 
of 1%. In addition, we choose the two coupling coeffi- 
cients g A /2-K = 472 MHz and 5 b /2tt = 549 MHz to dif- 
fer by approximately 15%. It is remarkable that, also 
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FIG. 5: (Color online) Robustness of quantum switch to fab- 
rication imperfections. Solid blue (dark grey) line: Numerical 
simulation of the quantum switch coupling coefficient |2(?!w | 
as a function of the frustration /,P C . Here, we assume a rel- 
atively large spread of 1% for the resonators center frequen- 
cies^ and a difference of approximately 15% between g A and 
g B . Dashed green (light grey) line: Plot of |2g!w| extracted 
from the generalized switch setting condition of Eq. (|18p for 
the same parameter spread as in the numerical simulations. 
For both the analytical and numerical result the switch set- 
ting condition is fulfilled (see black arrows). 

in this more general case, the switch setting condition 
can be fulfilled easily. We confirm this result by means 
of numerical simulations [solid blue (dark grey) line in 
Fig. [5] of the full Hamiltonian of Eq. (JTUJ), assuming fab- 
rication imperfections. Interestingly, in contrast to the 
case where A A = A B and g A — g B , we observe a dif- 
ferent behavior of the analytical and numerical curves of 
Fig. [5] when moving far away from the qubit degeneracy 
point. The reasons behind this fact rely on the conditions 
used to obtain the second-order Hamiltonian of Eq. (fT7|) . 

If ^AB ^ max {fAfB sin26l / 2A A'3A5B sin26l / 2A B:5AB}: 

as for the parameters chosen here, this Hamiltonian does 
not represent an accurate approximation anymore. In 
this expected, only a partial agreement between 

analytics and numerics is found. Nevertheless, a clear 
switch setting condition is obtained in both cases. We no- 
tice that the numerical switch setting condition is shifted 
towards the degeneracy point with respect to the analyt- 
ical solution. This is due to the detuning S AB present in 
Eq. (|17p . which is not accounted for when plotting the 
analytical solution. 

All the above considerations clearly show that the re- 
quirements on the sample fabrication are substantially 
relaxed. 



B. A quantum switch protocol 

We now propose a possible switching protocol based on 
three steps and discuss two different variants to shift from 
the zero-coupling to a finite-coupling condition charac- 



terized by a coupling coefficient g™. It is important to 
stress that this protocol is independent of the specific im- 
plementation (capacitive or inductive) of the switch. For 
definitcness, we choose a quantum switch based on a flux 
qubit in the following. 

(i) - First, we initialize the qubit in the groundstate |g). 

(ii) - Second, in order to fulfill the switch setting condi- 
tion, we choose the appropriate detuning A by changing 
the quasi-static bias of the qubit. For the switch oper- 
ation to be practical, we assume A = A 1 > 0. In this 
way, the sign of the coefficient in front of the <r 2 -operator 
of Eq. (|14p remains positive and, as a consequence, the 
switch is off in the groundstate |g), i.e., g'Jw = 0. 

(Hi) - Third, the state of the quantum switch can 
now be changed from off to on in two different ways, 
(a) or (b). 

(a) - Population-inversion. The qubit is maintained 
at the bias point preset in (ii). Its population is then 
inverted from |g) to |e), e.g., applying a Rabi 7r-pulse of 
duration t^. Such a pulse effectively changes the switch 

to the on-state, giw = 2g AB . In this case, g™ = 2g AB . 
Under these conditions, the two resonators are effectively 
coupled and the A-to-B transfer time is t = ir/2g°^ ll 
which also constitutes the required time-scale for most 
of the operations to be discussed in Subsec. IIII CI 

(b) - Adiabatic- shift pulse. We opportunely change the 
quasi-static bias of the qubit by applying an adiabatic- 
shift pulsed In this way, the qubit transition frequency 
becomes effectively modified. As a consequence, the de- 
tuning A is changed from A 1 to A 2 such that giw = 
g sw = g AB — g 2 sin 2 6/ A 2 ^ 0. In other words, the ge- 
ometric and dynamic coupling coefficients are not bal- 
anced against each other anymore and the switch is set 
to the on-state. In this case, = g sw . The rise 
time t risc of the shift pulse has to fulfill the conditional^ 
V<C > i ri se £ max {2tt/5 q , 

Variant (b) strongly benefits from the dependence of 
<7 SW on the external quasi-static bias of the qubit [see 
Figs. (3Jc) and[3Jd)]. We can distinguish between two 
possible regimes. The first regime is for a flux bias close 
to the qubit degeneracy point, where the qubit-resonator 
detuning is reduced and, thus, A 2 < A 1 . In this case, the 
dynamic contribution to g sw dominates over the geomet- 
ric one. This enables us to achieve very large resonator- 
resonator coupling strengths, which is a highly desir- 
able condition to perform fast quantum operations (e.g., 
cf. Sec. IIII C|) . The second regime is for a flux bias far 
away from the qubit degeneracy point, where the qubit- 
resonator detuning is increased and, thus, A 2 > A,. In 
this case, the geometric contribution to <7 SW dominates 
over the dynamic one. Since very far away from the qubit 
degeneracy point g sw — » |2<? AB | [cf. Subsec. IIII Al and 
Figs. QIa) anddjb)], operating the system in the second 
regime allows us to probe the pure geometric coupling 
between A and B. This would constitute a direct mea- 
surement of the geometric second-order coupling when 
m aq a W M Q Q » m. 
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Advanced applications: Nonclassical states 
and entanglement 



We now provide a few examples showing how the quan- 
tum switch architecture can be exploited to create non- 
classical states of the microwave radiation as well as en- 
tanglement of the resonators and qubit degrees of free- 
dom. In this subsection, when we discuss about the qubit 
we refer to the one used for the quantum switch opera- 
tion. If the presence of another qubit is required, we refer 
to it as the auxiliary qubit. 

Fock state transfer and entanglement between the res- 
onators. First, we assume the quantum switch to be 
turned off, e.g., following the protocol outlined in Sub- 
sec. IIII Bl with the qubit in the groundstate |g). In ad- 
dition, we assume resonator A to be initially prepared 
in a Fock state |1) A , while cavity B remains in the vac- 
uum state |0) B . Following the lines of Ref. 0, for ex- 
ample, a Fock state |1) A can be created in A by means 
of an auxiliary qubit coupled to it. A population inver- 
sion of the auxiliary qubit (via a 7r-pulse) and its subse- 
quent relaxation suffice to achieve this purpose. Then, 
we turn on the quantum switch for a certain time t fol- 
lowing cither one of the two variants (a) or (b) intro- 
duced in Subscc. HlIBI The initial states are |e) |1) A |0) B 
and |g) |1) A |0) B for (a) and (b), respectively. The quan- 
tum switch is now characterized by an effective coupling 
<7™ and the dynamics associated with the Hamiltonian of 
Eq. (Q3J) is activated for the time t. In this manner, a co- 
herent linear superposition of bipartite states containing 
a Fock state single photon^i 7 - ^ 60 ' 71 can be created 

cos (0) |1>A |0) B + e W2 CO) |0) A |1) B , (19) 

where the qubit state does not change and qubit and res- 
onators remain disentangled. If we choose to wait for a 
time t = 7r/2g™, we can exploit Eq. (|T9|) as a mecha- 
nism for the transferring of a Fock state from resonator 



A to resonator B, |1) A |0)j 



|0) A |1) B . In this case, 



also the resonators remain disentangled. It is notewor- 
thy to mention that the controlled transfer of a Fock 
state between two remote locations constitutes the basis 
of several quantum information devices^ 2 - If we choose to 
wait for a time t = ir/4g°£ instead, we can achieve max- 
imal entanglement between the two remote resonators. 
This goes beyond the results obtained in atomic systems, 
where two nondegenerate orthogonally polarized modes 
of the same cavity have been used to create mode entan- 
glement^ 

Tripartite entanglement and GHZ states. We follow a 
modified version of variant (a) of the switching proto- 
col. We start from the same initial conditions as in the 
previous example. Resonator A is in |1) A and resonator 
B is in |0) B . The qubit is in |g) and the switch setting 
condition is fulfilled, i.e., the switch is off. We then apply 
a 7r/2-pulsc to the qubit bringing it into the symmetric 
superposition 7 -^ (|g) + |c)) /v2- Then, the state of the 



system is still disentangled and can be written as 

|g)|l) A |0)B + |e)|l) A |0) B 
V2 



(20) 



Now, the Hamiltonian of Eq. (Q3J) yields the time evolu- 
tion 



i(|g)|l) A |0) B + cos(0)|e)|l) A |0) B 
+ e^/ 2 sin(0)|e) |0) A |1) B ) 



(21) 



for the state of the quantum switch. Under these con- 
ditions, the dynamics displayed in Eq. (f2~Tj) is character- 
ized by two distinct processes. The first one acts on the 

|g> |!)a I°)b P art of the initial state 01 E q- ©■ This 
process is actually frozen because the quantum switch is 
turned off when the qubit is in |g). The second process, 
instead, acts on the |e) |1) A |0) B part of the initial state, 
starting the transfer of a single photon from resonator A 
to resonator B and vice versa. If during such evolution 
we wait for a time t = 7r/2g™, a tripartite entangled 
state 



|g>|l> A |0> B +e W2 |e) |0) A |1) B 



V2 



(22) 



of the GHZ classic is generated. Here, the two res- 
onators can be interpreted as photonic qubits because 
only the Fock states |0) AB and |1) AB are involved. 
Hence, Eq. (f2"2"| represents a state containing maximal 
entanglement for a three-qubit system, which consists 
of two photonic qubits and one superconducting (charge 
or flux) qubit. The generation of GHZ states is impor- 
tant for the study of the properties of genuine multi- 
partite entanglement. Interestingly, the quantum na- 
ture of our switch is embodied in the linear superposi- 
tion of |g) |1) A |0) B and |e) |1) A |0) B of the initial state of 
Eq. CD). 

Entanglement of coherent states. Finally, we show how 
to produce entangled coherent states of the intracavity 
microwave fields of the two resonators. These are pro- 
totypical examples of the vast class of states referred to 
as Schrodingcr cat states i 21 i 43 ' 75 i 76 ' 77 This time, we start 
with cavity A populated by a coherent state |a) A instead 
of a Fock state |1) A . Again, cavity B is in the vacuum 
state |0) B and the qubit in the symmetric superposition 
state (|g) + |e)) /V2, i.e., a modified version of variant 
(a) of the switching protocol is again employed. The 
total disentangled initial state can be written as 



5> |Q) A |0) b + |e> 10), 
V2 



(23) 



The resulting dynamics associated with the Hamiltonian 
of Eq. (fT4)) yields a time evolution similar to that shown 
for Fock states in Eq. (f2"T]) . In this case, the part of the 
evolution involving |c) can be calculated either quantum- 
mechanically or employing a semi-classical model. In 
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both cases, after a waiting time t = ir/2g™, resonator 
B is in the state \a) B and A in the vacuum state |0) A . 
However, Eq. (f23|) contains an initial linear superposition 
of |g) |a) A |0) B and |e) |a) A |0) B , requiring a quantum- 
mechanical treatment of the time evolution. From this, 
one finds that, after the waiting time t = ir/2g™ the 
quantum switch operation creates the tripartite GHZ- 
type entangled state 



10)! 



o3V 



a 



V2 



(24) 



where ip is an arbitrary phase. Again, the creation of such 
states clearly reveals the quantum nature of our switch, 
showing a departure from standard classical switches^ 
Remarkably, the state of Eq. ([24]) describes the entangle- 
ment of coherent ("classical") states in both resonators. 
This feature is peculiar to our quantum switch and can- 
not easily be reproduced in atomic systems^ In prin- 
ciple, in absence of dissipation the quantum switch dy- 
namics continues transferring back the coherent state to 
cavity A. In order to stop this evolution, an ultimate 
measurement of the qubit along the x-axis of the Bloch 
spher o 7 ' 78 is necessary. This corresponds to a projection 
associated with the Pauli operator a x , which creates the 
two-resonator entangled state 



|a) A |0)j 



o3<P 



10), 



V2 



(25) 



This state is decoupled from the qubit degree of freedom. 

Obviously, all the protocols discussed above need suit- 
able measurement schemes to be implemented in reality. 
For instance, it is desirable to measure the transmitted 
microwave field through both resonators and, eventually, 
opportune cross-correlations between them by means of a 
double heterodyne detection apparatus similar to that of 
Ref. l60l . In addition, a direct measurement of the qubit 
state, e.g., by means of a DC SQUID coupled to it&i2 
would allow for the full characterization of the quantum 
switch device. 

In summary, we show that a rich landscape of nonclas- 
sical and multipartite entangled states can be created 
and measured by means of our quantum switch in two- 
resonator circuit QED. 



IV. TREATMENT OF DECOHERENCE 

The discussion in the previous sections implicitly as- 
sumes pure quantum states. In reality, however, a quan- 
tum system gradually decays into an incoherent mixed 
state during its time evolution. This process, called de- 
coherence, is due to the entanglement of the system with 
its environment and it is known to be a critical issue for 
solid-state quantum circuits. Since it is difficult to decou- 
ple these circuits from the large number of environmental 
degrees of freedom to which they are exposed, their typ- 
ical decoherence rates cannot be easily minimized. Usu- 



ally, they are in the range from 1 MHz to 1 GHz, depend- 
ing on the specific implementation. In this section, we 
discuss the impact of the three most relevant decoherence 
mechanisms on the quantum switch architecture. These 
are: First, the population decay of resonators A and B 
with rates k a and k b , respectively; second, the qubit 
relaxation from the energy excited state to the ground- 
state at a rate 7 r due to high-frequency noise; third, the 
qubit dephasing (loss of phase coherence) at a pure de- 
phasing rate 7^ due to low-frequency noise. We show 
by means of detailed analytical derivations that, despite 
decoherence mechanisms, a working quantum switch can 
be realized with readily available superconducting qubits 
and resonators. 

Decoherence processes are most naturally described in 
the qubit energy eigenbasis {|g) , |e)}. Under the Markov 
approximation, the time evolution of the density matrix 
of the quantum switch Hamiltonian H of Eq. pT|) is de- 
scribed by the master equation 



pH) 



(26) 



n=l 



Here, C n is the Lindblad superoperator defined as L n p = 

7 „ (X„plt - XlX n p/2 - pXlXji) . The indices n = 

1,2,3,4 label the decay of resonator A, the decay of 
resonator B, qubit relaxation, and qubit dephasing, re- 
spectively. Consequently, the generating operators are 



J = At 



X, 



b,xl = tf,x a = a-,x. 



n 



<7 Z . The corresponding decoherence 

K A> 72 = «B) 73 = Tr, and 74 = %/ 2 - For 



X l = a, X\ 
and A 4 = 
rates are 7 1 

the resonators, n A and k b are often expressed in terms 
of the corresponding loaded quality factors Q A and Q B , 
k a = lo a /Q a and n B = uj b /Q b , respectively. Although 
in general all four processes coexist, in most experimen- 
tal situations one of them dominates over the others. In 
fact, it is a common experimental scenario that 7^ ^> 7 r , 
for example in the special case of a flux qubit operated 
away from the degeneracy point (see, e.g., Ref. I58T ). In 
this situation, we can extract pessimistic relaxation and 
dephasing rates from the literature ; 58 ' 79 ' 80 ' 81 7 r — 1 MHz 
and 7^ ~ 200 MHz. In other words, dephasing is the 
dominating source of qubit decoherence^ The decay 
rates of the resonators can be engineered such that 7 - - 
K A; k b ~ 7r ^ Tip- For these reasons, hereafter we focus 
on dephasing mechanisms only. Hence, we analyze the 
following simplified master equation 



1 



Hp - pH) 



(27) 



where C^p = C A p = {-/ ip /2)(aja z - p). 

The impact of qubit dephasing on the switch operation 
depends on the chosen protocol (cf. Subsec. HlIB[) . When 
employing the population-inversion protocol, qubit de- 
phasing occurs within the duration time t n of the control 
7r-pulses. The time t ff coincides with the inverse of the 
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qubit Rabi frequency and can be reduced to less than 
1 ns using high driving power In this way, the time 
window during which the qubit is sensitive to dephasing 
is substantially shortened. However, it is more favorable 
to resort to the adiabatic-shift pulse protocol. In this 
case, the qubit always remains in |g) resulting in a com- 
plete elimination of pulse-induced dephasing. The rele- 
vant time scale during which dephasing occurs is there- 
fore set by the operation time of the switch between two 
on-off events. Naturally, this time should be as long as 
possible if we want to perform many operations. 

The effect of dephasing during the switch operation 
time is better understood by inspecting the effective 
quantum switch Hamiltonian H eS of Eq. (fT4|) . In Sub- 
scc. UlIAl we deduce this effective Hamiltonian by means 
of a Dyson series expansion. This approach is very 
powerful and compact when dealing with the analy- 
sis of a unitary evolution. However, when treating 
master equations, wc prefer to utilize a variant of the 
well-known Schricffer- Wolff unitary transformation ; 2 ^ 84 
UHU\ where 



U = exp 



— 1 " a' — a a 



A (* 



(28) 



and U* is its Hcrmitian conjugate. In the large-detuning 
regime, g A sin 9, c/ B sin 9 <C A, we can neglect all terms 
of orders (g A sin#/A) 2 , (g B sin#/A) 2 , g A g B sin 2 0/A 2 , or 
higher. After a transformation into an interaction picture 
with respect to the qubit and both resonators (cf. Sub- 
scc. MI A[) and performing opportune RWAs, we obtain 
again H cS of Eq. (fT?)) . The master equation govern- 
ing the time evolution of the effective density matrix 



P 



eflf 



U pU< then becomes 



ieff 



^H cS ) +Cfp cS . 



(29) 



The analysis is complicated by the fact that also the Lind- 
blad superoperator Cp has to be transformed. For the 
sake of simplicity, we can assume g A = g B = g and find 



Soff -off 



g sin ( 



(30) 



When deriving this expression, all terms of O (g sin 9/ A) 
are explicitly neglected by a RWA. This approximation 
relies on the condition (7 /A) gsin^ <C A, which is well 
satisfied in the large-detuning regime as long as 7 < A. 
The latter requirement can easily be met by most types 
of existing superconducting qubits. In the frame of H cS , 
C v p eS has the standard Lindblad dephasing structure 
and the qubit appears only in the form of o^-operators. 
Since the initial state of the switch operation is charac- 
terized by cither no (adiabatic-shift pulse protocol) or 



only very small (population-inversion protocol) qubit co- 
herences, the effect of Cpf) cS on the time evolution of 
the system is negligible. All other nonvanishing terms 
are comprised in the expression 27^ x 0[(g sintf/A) 2 ] of 
Eq. (f3"0"l) and scale with a factor smaller than 7° ff = 

^■lip {9 sm 9/ A) 2 . Hence, the operation of the quantum 
switch is robust to qubit dephasing on a characteristic 
time scale T^ ff = l/7^, ff l/7„. For completeness, it 
is important to mention that the higher-order terms of 
Eq. (|30[) can contain combinations of operators such as 
a^a and wb. In this case, T° would be reduced for a 
large number of photons populating the resonators. For- 
tunately, this does not constitute a major issue since the 
most interesting applications of a quantum switch require 
that the number of photons in the resonators is of the or- 
der of one. 

In summary, we show that for suitably engineered cavi- 
ties the quantum switch operation time for the adiabatic- 
shift pulse protocol is limited only by the effective qubit 
dephasing time T° ff . The latter is strongly enhanced with 
respect to the intrinsic dephasing time T v = 1/7^. In 
this sense, the quantum switch is superior to the dual 
setup, where two qubits are dispersively coupled via one 
cavity bus£ Moreover, the intrinsic dephasing time T 
and, consequently, T° are further enhanced by choos- 
ing a shift pulse which sets the on-state bias near the 
qubit degeneracy point.— As explained in Subsec. MI Bl 
this regime takes place for a qubit-resonator detuning 
A 2 < A x . In this case, the switch coupling coefficient 
is also substantially increased because of a dominating 
dynamic interaction. As a consequence, this option is 
particularly appealing in the context of the advanced ap- 
plications discussed in Subscc. MI CI Finally, wc notice 
that for the population-inversion protocol the switch op- 
eration time could be limited by the qubit relaxation time 
T r = l/7 r However, the switch setting condition is typ- 
ically fulfilled for a bias away from the qubit degeneracy 
point. There, T r is considerably enhanced by both a re- 
duced^ 8 - sin 9 and by the Purcell effect of the cavities.— 



V. AN EXAMPLE OF TWO-RESONATOR 
CIRCUIT QED WITH A FLUX QUBIT 

In this section, we focus on the geometry sketched 
in Fig. QJc) and present one specific implementation of 
two-resonator circuit QED. As a particular case, the de- 
scribed setup can be operated as a superconducting quan- 
tum switch. In this example, we consider microstrip res- 
onators. Coplanar wave guide resonators can also be used 
without significantly affecting our main results. In addi- 
tion, we choose a flux qubit because this is our main topic 
of research] 12 ' 58 i 60 ' 81 Moreover, as shown in Subsec. MI A[ 
the dynamic properties of the quantum switch are inde- 
pendent of specific implementations. As a consequence, 
in this section we concentrate on its geometric properties 
only. It is worth mentioning again that such properties 
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FIG. 6: (Color online) A possible setup for two-resonator circuit QED with a flux qubit. (a) Overall structure (dimensions 
not in scale). Two microstrip resonators A and B (thick blue lines) of length £ m simultaneously coupled to a flux qubit loop 
[magenta (middle grey) rectangle]. C a in , C b in , C a out , and C b out : Input and output capacitors for A and B. The dashed black 
box indicates the region of the close-up shown in (b). £ a \ m - Length of the region used for the FASTHENRY— simulations, 
(b) Close-up of the region which contains the flux qubit loop in (a). ^ q i and £ q 2'- Qubit loop lateral dimensions. W^. Width 
of the qubit lines. d mq : Distance between the qubit and each resonator. The dashed black line denominated as S marks the 
cross-section reported on the top part of the panel. t q : Thickness of the qubit loop lines, (c) As in (b), but without the qubit 
loop. W m and t m : Width and thickness (see cross-section S) of the two microstrip resonators. H s : Height of the dielectric 
substrate. The reference axis Oz is also indicated (cf. Appendix [B]). Both in (b) and (c), a in , a out , b in , and b out represents the 
input and output probing ports used in the simulations, (d) Current density distribution at high frequency (5 GHz) for the 
structures drawn in (c). The currents are represented by small arrows, green (light grey) for resonator A and blue (dark grey) 
for resonator B. (e) Current density distribution at high frequency (5 GHz) for the structures drawn in (b). The two black 
arrows indicate two high-current-density channels between the two resonators. The dashed black box marks the close-up area 
shown in (f). (f) Close-up of one of the two geometric second-order interaction channels. 



are inherent to circuit QED architectures and constitute 
a fundamental departure from quantum optical systems. 

In Figs. Uta) and[6jb), the design of a possible two- 
resonator circuit QED setup is shown. The overall 
structure is composed of two superconducting microstrip 
transmission lines, which are bounded by input and out- 
put capacitors, C a in , C b in , C a out , and C h out . This ge- 
ometrical configuration forms the two resonators A and 
B. The input and output capacitors also determine the 
loaded or external quality factors Q A and Q B of the two 
resonators.— Both A and B are characterized by a length 
£ m , which defines their center frequencies / A and / B . We 
choose £ m ~ A m /2 ~ 12 mm, where A m = A A = A B 
is the full wavelength of the standing waves on the res- 



onators. The superconducting loop of the flux qubit cir- 
cuit is positioned at the current antinode of the two A m /2 
resonators. 

In Fig. |B(c), only the two microstrip resonators A and 
B are considered. They are chosen to have a width 
W m = 10 urn and a thickness t m = 100 nm. The height of 
the dielectric substrate between each microstrip and the 
corresponding groundplane is H s = 12.3 urn. The sub- 
strate can opportunely be made of different materials, 
for example silicon, sapphire, amorphous hydrogenated 
silicon, or silicon oxide, depending on the experimen- 
tal necessities. A detailed study on the properties of 
a variety of dielectrics and on the dissipation processes 
of superconducting on-chip resonators can be found in 
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Refs. mi3liilli and [H The aspect ratio W m /H s is 
engineered to guarantee a line characteristic impedance 
Z c ~ 50 £1, even if this is not a strict requirement for the 
resonators to function properly^ The remaining dimen- 
sions of our system are shown in Fig. [6^b): The lateral 
dimensions t ql = 200 |xm and £ q2 = 87 |J.m of the qubit 
loop, the width W q = 1 urn of each line forming the qubit 
loop, the interspace d mq = 1 |^m between qubit and res- 
onators, and the thickness i q (= t m ) = 100 nm of the 
qubit lines. The dimensions of the qubit loop are cho- 
sen to optimize the qubit-resonator coupling strengths. 
This geometry results in a relatively large inductance 
L q ~ 780 pH (number obtained from FASTHENRY sim- 
ulations^ cf. Table |TJ|. Despite the large self- inductance 
L q , reasonable qubit coherence times are expected (see, 
e.g., Refs. [H and HI). Moreover, in the light of Sec. HVl 
these coherence times easily suffice for the operation of a 
superconducting quantum switch, where the qubit acts as 
a mere mediator for the exchange of virtual excitations. 

In our numerical simulations (cf. Appendix [Bj) , we 
restrict ourselves to the region indicated by the black 
dashed box in Fig. Ela), the close-up of which is shown 
in Fig. HJb) and, in the absence of the flux qubit loop, 
in Fig. Etc). This region is characterized by a length 
£ sim = 500 |a.m of the resonators and is centered at a posi- 
tion where the magnetic field reaches a maximum (antin- 
ode) and the electric field reaches a minimum (node). We 
notice that magnetic and electric fields can equivalcntly 
be expressed in terms of currents and voltages, respec- 
tively. There are two main hypotheses behind the valid- 
ity of our simulation results for the entire two-resonator- 
qubit system. These are the uniformity of the electro- 
magnetic field (voltage and current) in the simulated 
region and the abruptly 9 ^ increasing geometric distance 
between resonators A and B outside of it [see sketch of 
Fig. Ela)]. The three main implications of the above as- 
sumptions are explained in the following. First, all cou- 
pling strengths are dominated by inductive interactions 
and there are no appreciable capacitive ones. Inside the 
simulated region, in fact, the voltage is practically char- 
acterized by a node, which results in a vanishing coupling 
coefficient. Outside the simulated region, the effective 
distance d eS between the cavities strongly increases to- 
gether with the geometric onei 95 ' 96 As a consequence, the 
geometric first-order capacitance coal/ d eff becomes neg- 
ligible. Second, the coupling coefficients between qubit 
and resonators can be obtained without integrating over 
the spatial distribution of the mode. This is because 
of the uniformity of the field, which, for all practical 
purposes, is constant over the restricted simulated re- 
gion. Third, the geometric first-order coupling between 
the two resonators, which is proportional to their mutual 
inductance m, is still accurately determined. In fact, 
outside the simulated region any additional contribution 
to m becomes negligible. For all the reasons mentioned 
above, we are allowed to use the FASTHENRY 9 ^ cal- 
culation software for our simulations. In this section, 
we utilize two different versions of FASTHENRY, one 



for superconducting materials and one for almost perfect 
conducting ones. We use the second version only when 
we want to obtain current density distributions or the 
frequency dependence of an inductance. In these cases, 
due to technical limitations of the software, we cannot 
use the version valid for superconductors? 9 ^ 

Figures ^d) and [6je) display the simulated cur- 
rent density distributions at a probing frequency of 
5 GHz (high-frequency regime) for the different struc- 
tures drawn in Fig. [He) and [DJb) , respectively. Simi- 
lar results can be found in a range between 1 GHz and 
10 GHz (data not shown). Without loss of generality, 
these simulations arc performed for almost perfect con- 
ductors using a FASTHENRY version which does not 
support superconductivity. The results of Fig. E^d) 
clearly show that the two microstrip lines are regions 
characterized by a high current density separated by a 
region with a low current density in absence of the flux 
qubit loop. In this case, the geometric interaction be- 
tween resonators A and B is reduced to a bare first-order 
coupling, which turns out to be very weak. On the con- 
trary, in Fig. [He) the presence of the qubit loop clearly 
opens two new current channels between A and B. These 
are located at the position of the upper and lower qubit 
loop segments of Fig. EJb). For clarity, the close-up of 
one of these channels is shown in Fig.[6]Jf). Notably, the 
two additional current channels of Fig.^e) represent the 
geometric second-order coupling. 

We now study in more detail the relationship between 
geometric first- and second-order inductances for the 
structures of Fig. (SJb) and^c). The notation is that of 
Subsec. fll CI and Figs. [2jd)-(f). All quantities are com- 
puted numerically with the aid of FASTHENRY for su- 
perconducting materials^ assuming a London penetra- 
tion depth A L = 180 nm. We notice that, in this case, 
the simulated inductances are independent of the prob- 
ing frequency. In a first run of simulations, we calcu- 
late pure first-order inductances only (cf. Appendix IB"]) . 
These are the simulated test inductances L* a and L* h 
from which we obtain the self-inductances L ra and L rh of 
resonators A and B (more details in the next paragraph), 
the first-order mutual inductance m between the two res- 
onators, the self-inductance L q of the qubit loop, and the 
mutual inductances M qa and M qh between qubit and res- 
onators. In a second run of simulations, we calculate di- 
rectly (cf. Appendix iBj) the sum of first- and second-order 
inductances. These are the renormalized test inductances 
L* a and L* h of the portions of resonators A and B shown 
in Fig. Etb) and the total mutual inductance fh between 
the two resonators. The difference fh — rn = 4.78192 pH, 
i.e., the geometric second-order coupling, coincides up to 
the sixth significant digit with the quantity M qa M qh / L q 
expected from our general three-node network approach 
of Subsec. IIIBI and, equivalcntly, from the three-circuit 
theory of Subsec. HTC] (cf. Table [Q). We also find that 
the dominating geometric coupling between A and B is 
not the first-order inductance m = 2.90130 pH, but the 
second-order one. The ratio between second-order and 
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TABLE I: Relevant parameters for a possible two-resonator circuit QED setup based on a superconducting flux qubit. The 
various constants are described in Subsec. [V] and reported in Fig. [2] All inductances are simulated using the version of FAS- 
THENRY for superconducting materials.— The capacitances are calculated analytically. All geometric second-order inductances 
are computed analytically and numerically and then compared to each other for consistency. We find an excellent agreement 
between the estimates M qa M qb /L q and fh — m for the second-order mutual inductance. Also, the shift inductances M qa /L q 
and M qb /L q coincide with their counterparts L sa = L* a — Z* a and L sb = L* b — Z* b , respectively. These parameters are suitable 
for the implementation of a superconducting quantum switch. 
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first-order inductances is (m — m)/m ~ 1.6. In addition, 
the numerical simulations yield the two shift test induc- 
tances |Z* a - L*J = \L* h - L* h \ = 4.78100pH (cf. also 
next paragraph). These shifts renormalize the bare cen- 
ter frequencies / A and / B of resonators A and B, respec- 
tively, and are found to be in very good agreement up 
to several decimal digits with their analytical estimates 
L sa = M qa /L q and L sh = M qh /L q of our three-circuit 
theory of Subsec. Hi Cl fcf. Table|J). We point out that, in 
our definition, the quantities L aa and L sb are strictly pos- 
itive. Remarkably, our simulations reveal that L* a > L*. a 
and L* h > L* b , reproducing the minus sign in the ex- 
pressions L ra — L sa and L rb — L sh of Fig. HJf). These 
findings confirm the superiority of the three-circuit the- 
ory of Subsec. [II Cl over the simple model which results in 
the Hamiltonian of Eq. ([9]). In the case of purely induc- 
tive interactions, the resonators suffer a small blue-shift 
of their center frequencies, i.e., a shift towards higher val- 
ues. This is opposite to the redshift, i.e., towards lower 
frequencies, experienced by the resonators for a pure ca- 
pacitive coupling [cf. Subsec. [II CI and Fig. [2(c)]. 

The numerical values of all parameters discussed above 
are listed in Table fl] The values of the bare self- 
inductances of the two resonators are first evaluated for 
the test length l sim in absence of the qubit loop. This 
yields the simulated test inductances L* a and L* h . Then, 
L* a and L* h are extrapolated to the full length £ m of 
each microstrip resonator to obtain L ra and L rb , respec- 
tively. In the presence of the qubit loop, the simulated 
test inductances L* a and L* h can also be found. The 
resonator capacitances per unit length, c ra and c rb , are 
calculated analytically by means of a conformal mapping 
technique^ 

/ 8H W cfi \ 
c ra = c rh = 2ne e r ln^ + ^-J . (31) 

Here, e = 8.854 x 10~ 12 F/m is the permittivity of 



vacuum (electric constant)^ e r = 11.5 the relative di- 
electric constant of the substrate (in our example, sap- 
phire or silicon; other dielectrics could be used), and 
W£ ff = W m + 0.398t m [l + ln(2HJt m )] the effective width 
of the resonators i2£ As a consequence, the resonator ca- 
pacitances are C ra = £ m c ra and C rb = i m c xh . Finally, 
from the knowledge of the velocity of the electromagnetic 
waves inside the two resonators, c A = i m / \/ (L la C ra ) and 
c b = £ m / y/ (i rb C rb ), one can find the full wavelengths 
A A = c A // A and A B = c B // B of the two resonators. As 
before, all these results are summarized in Table HI 

We now analyze the frequency dependence of the ge- 
ometric first- and second-order coupling coefficients, i.e., 
the first- and second-oder mutual inductances, for a 
broad frequency span between 1 Hz and 10 GHz. Again, 
we assume almost perfectly conducting structures and 
use the FASTHENRY version which does not support su- 
perconductivity. The results are plotted in Figs. [7] and [8] 
In Fig.[7fa), we plot the frequency dependence of the sim- 
ulated inductances L q (which is renormalized by a factor 
of 8.5 for clarity) and M qa = M qb . From these, we then 
compute the expression M qa M qh / L q for the second-order 
mutual inductance as derived in Subsecs. Ill Bl and III CI 
This expression is plotted in Fig. [TJb). In Fig. \T[c), we 
plot the bare second-order mutual inductance fh — m as a 
function of frequency. Remarkably, comparing Fig. [Tfb) 
to Fig. 17(c) , we find M qa M qh / L q = fh — m with very 
high accuracy over the entire frequency range. In the 
frequency region of interest for the operation of a quan- 
tum switch, i.e., from approximately 1 GHz to 6 GHz, we 
find L q ~ 63.02 pH, M qa = M qh ~ 7.37 pH, and, con- 
sequently, M qa M qh /L q = fh — m ~ 7.33 pH. All these 
results prove again the general validity of the derivations 
of Subsecs. [HB] and HTCl 

Finally, we study the scattering matrix elements be- 
tween resonators A and B both without and with flux 
qubit loop. In absolute value, these elements correspond 



19 



64.0 



£63.5 



CD 
O 
£= 

cd 

o 63.0 
c 



62.5 - 



a L /8.5 

v ; q 

M = M J 

qa v qb' 



10" 



10 3 

frequency (Hz) 



10' 



x 



7.3- 



- ra 7.2- 



(b) 




<E 7. 



10 3 

frequency (Hz) 



10' 



FIG. 7: (Color online) FASTHENRY simulation results for the frequency dependence of some relevant first- and second-order 
inductances relative to our example of two-resonator cicruit QED. Vertical dashed black lines: Frequency region of interest for 
the operation of the quantum switch from 1 GHz to 6 GHz. (a) Magenta (middle grey) line: Qubit loop self-inductance L q 
renormalized by a factor of 8.5 for clarity. Green (light grey) line: Mutual inductance M qa ( = M qb ) between the qubit and 
resonator A (or B). (b) Bare second-order mutual inductance between the two resonators calculated with the results from (a) 
using the expression M qa M qb /L q . (c) Bare second-order mutual inductance between the two resonators fh — m. The agreement 
with (b) is excellent. 



to the isolation coefficients between A and B. As before, 
the FASTHENRY simulations are performed within the 
regions of Figs.[SJc) and[HJb), respectively. In these fig- 
ures, we also define the input and output probing ports 
used in the simulations as a in , a out , b in , and b out , respec- 
tively. Under these assumptions, the scattering matrix 
element S ab = S ba in absence of the flux qubit loop is 
given b y 95 i 96 



5 ab ^201og 



j+=o 



m 

= 20 log — 



(32) 



where I b is a test current wave incident on the input 
probing port b in of resonator B. The current — 1~ cor- 
responds to the outgoing wave from the input probing 
port a in of resonator A. The remaining current waves 
incident on the ports of the two resonators are / + = 

1 1£ , I a t , I b J. In a similar way, the scattering ma- 
trix element S ab = S ba in presence of the flux qubit loop 
is given by 



S. 



ab 



20 log 



m 



(33) 



We note that the same results as in Eqs. (|32|) and (|33|) 
are obtained replacing the input probing port b in with 
the output probing port b out for the incident wave. In 
this case, the associated current I b has to be exchanged 

with I b . Similar substitutions apply for the probing 
port and associated current of the outgoing waves. In the 
literature^ the outgoing waves are often denominated as 



reflected waves. Equation (|32"|) can be straightforwardly 
found via the definitions of mutual and self-inductance, 
ml b = <I> ba = L* a I~ . There, $ ba is the flux generated 
in the portion of resonator A by the current flowing in 
the portion of resonator B of Fig. EJc). Similar argu- 
ments lead to Eq. (f3"3)) . When considering superconduct- 
ing materials, the scattering matrix elements between A 
and B without and with flux qubit loop can be evalu- 
ated inserting the opportune numbers reported in Table U 

- 38.80 dB 



ab 



into Eqs. (32]) and (33]). This yields S. 
and S ab ~ — 30.18dB. If we want to calculate the scat- 
tering matrix elements between A and B over a broad 
frequency span (e.g., from 1 Hz to 10 GHz), we can con- 
sider again almost perfect conducting structures. In this 
case, the results are plotted in Figs. E[a) and[S]Jb). In 
the high frequency region from 1 GHz to 6 GHz, we find 
5 ab ~ - 37.66 dB and S ah ~ - 30.54 dB. These num- 
bers are in good agreement with the results obtained for 
superconducting materials. In addition, it is noteworthy 
to mention that the scattering matrix elements between 
A and B calculated here with FASTHENRY for almost 
perfect conducting structures are in excellent agreement 
with those evaluated for similar structures by means of 
more advanced software based on the method of mo- 



ments 



In conclusion, we study a detailed setup of two- 
resonator circuit QED based on a superconducting flux 
qubit. In this case, we prove that the geometric second- 
order inductance found with our three-node network ap- 
proach agrees well with that obtained from numerical 
simulations. Moreover, we give a set of parameters (many 
sets can easily be found) for which the second-order in- 
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FIG. 8: (Color online) FASTHENRY simulation results for 
the frequency dependence of the scattering matrix elements 
between resonators A and B considering almost perfect con- 
ducting structures. Frequency span: From 1 Hz to 10 GHz. 
(a) Scattering matrix element S ab in absence of the flux qubit 
loop, (b) Scattering matrix element S ab in presence of the 
flux qubit loop. Owing to the significant second-order mutual 
inductance between A and B, we find |S ab | < |5 ab |. 



ductancc dominates over the first-order one. 



We subsequently demonstrate the possibility of balanc- 
ing a geometric coupling against a dynamic second-order 
one. In this way, the effective interaction between the 
two resonators can be controlled by means of an exter- 
nal bias. Based on this mechanism, we propose possible 
protocols for the implementation of a quantum switch 
and outline other advanced applications, which exploit 
the presence of the qubit. 

Remarkably, we find that the quantum switch opera- 
tion is robust to decoherence processes. In fact, we show 
that the qubit acts as a mere mediator of virtual exci- 
tations between the two resonators, a condition which 
substantially relaxes the requirements on the qubit per- 
formances. 

Finally, we give detailed parameters for a specific setup 
of two-resonator circuit QED based on a superconducting 
flux qubit. We perform numerical simulations of the ge- 
ometric coupling coefficients and find an excellent agree- 
ment with our analytical predictions. In particular, we 
confirm the existence of a regime where the geometric 
second-order coupling dominates over the first-order one. 

In conclusion, our findings show that, in circuit QED, 
the circuit properties of the system are crucial to pro- 
vide a correct picture of the problem and also constitute 
a major difference with respect to atomic systems. This 
peculiar aspect of ciruit QED makes it a very rich envi- 
ronment for the prediction and experimental implemen- 
tation of unprecedented phenomena. 



VI. SUMMARY AND CONCLUSIONS 

In this work, we first introduce the formalism of two- 
resonator circuit QED, i.e., the interaction between two 
on-chip microwave cavities and a superconducting qubit 
circuit. Starting from the Hamiltonian of a generic three- 
node network, we show that the qubit circuit mediates 
a geometric second-order coupling between the two res- 
onators. For the case of strong qubit-resonator coupling, 
the geometric second-order interaction is a fundamen- 
tal property of the system. In contrast to the geomet- 
ric first-order coupling between the two resonators, the 
second-order one cannot be arbitrarily reduced by means 
of proper engineering. 

With the aid of two prototypical examples, we then 
highlight the important role played by circuit topology in 
two- resonator circuit QED. Our models reveal a clear de- 
parture from a less detailed theory based on the Hamilto- 
nian of a charge quantum circuit (e.g., a Cooper-pair box 
or a transmon) or a flux quantum circuit (e.g., a super- 
conducting one- or three- Josephson-j unction loop) cou- 
pled to multiple quantized microwave fields. We demon- 
strate that this simplified approach easily produces ar- 
tifacts. We also show that our three-node network ap- 
proach suffices to obtain correct results when including 
topological details appropriately into the definitions of 
the nodes. 
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APPENDIX A: HIGHER-ORDER CORRECTIONS 
TO THE CAPACITANCE AND INDUCTANCE 
MATRICES 



In Subscc. Ill B[ we account for corrections up to 
second-order capacitive and inductive interactions be- 
tween the elements of a three-node network. Throughout 
this work, we show that for a three-node network the ge- 
ometric second-order coupling coefficients can dominate 
over the first-order ones. For this reason, in the follow- 
ing we can safely assume vanishing first-order coupling 
coefficients, c = m = 0. Nevertheless, we notice that our 
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results would not be qualitatively affected even in the 
presence of appreciable first-order couplings. 

In this appendix, we demonstrate that third- and 
fourth-order capacitances and inductances are negligible. 
We start with the case of third-order corrections. There 
are two possible excitation pathways giving rise to third- 
order coupling coefficients. These pathways are between 
resonator A and qubit Q, A-Q-B-Q, or between resonator 
B and qubit Q, B-Q-A-Q. Assuming the two resonators 
to have identical properties, we only study the A-Q-B-Q 
pathway. In this case, from the knowledge of the electro- 
magnetic energy we can derive 
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where the inverse paths Q-A and Q-B-Q-A are already 
included. In the equation above, resonator B is only vir- 
tually excited. In the same equation, we identify the 
capacitance and inductance matrix elements up to third 
order 
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In circuit theory, it is well-known that the squares 
of the electromagnetic coupling coefficients are^ 
Cqb/^qqCbb < 1 and Af2 B /M QQ M B? < 1. This 
implies that the pure third-order capacitance and in- 
ductance are always smaller than the first-order ones, 



C ( aq ~ C AQ < C AQ and M$ - M Aq < M AQ . For typical 
experimental parameters, we find third-order processes to 
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be negligible, C AQ - ^ AQ ^ ^ AQ ^ ,„ AQ 
M A q. For example, using the parameters given in Scc.lVl 
yields M% B /M QQ M BB ~ 7.88 x 10~ 4 < 1. 

In a similar way, the fourth-order coupling coefficients 
for the excitation pathways A-Q-B-Q-B and, equiva- 
lently, B-Q-A-Q-A can easily be found. In this case, it 
is the qubit to be only virtually excited. The capaci- 
tance and inductance matrix elements up to fourth order 
become 
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When comparing the above equations to Eqs. (|A2[) and 
(|A3|) , respectively, we find that fourth-order processes are 
negligible for typical experimental parameters. 

In the light of all these considerations, all higher-order 
coupling coefficients can safely be ignored within the 
scope of this work. 



APPENDIX B: 



DETAILS OF THE FASTHENRY 
SIMULATIONS 



In this appendix, we discuss the details of our FAS- 
THENRY simulations^ First, we verify our hypothesis 
on the uniformity of the AC currents (corresponding to 
the magnetic fields) flowing on the resonators in the re- 
gions of Figs. HJb) and[5Jc). To this end, we derive the 
quantized current on one of the two resonators (e.g ., A) 
following similar calculations as in Refs. [22l and l60l 



4a OM) = i 



AO 



cos 



(Bl) 



where z represents a coordinate along the longitudinal di- 
rection of the resonator [see Fig. [DJc)] and t is the time. 
In Eq. (|Bip . the bosonic field operators are expressed in 
the Heisenberg picture. We notice that Eq. (|Bip is valid 
for the first mode of the A/2-resonator(s) considered in 
our example. The contribution from the second mode is 
negligible for two main reasons. First, the current is char- 
acterized by a node at the flux qubit loop position cho- 
sen here. Second, the qubit-rcsonator detuning becomes 
substantially larger, hence resulting in a correspondingly 
reduced coupling. The contribution form higher modes 
can also be neglected because of the increasing detuning. 

Substituting the numbers of Table |l] into Eq. (jBlj) and 
setting z = =F£ sim /2, we find that the two currents at 
the boundaries =F£ sim /2 are about 0.998i AO , where z A0 
is the maximum amplitude of the quantized current in 
Eq. (jBlj) . This maximum is obtained at the position 
z = of the mode antinode. The main implications 
of current uniformity over £ sim are explained in detail 
in Sec. |V1 In a similar way, we can also estimate the 
voltage contribution for the first mode at the boundaries 
=F£ sim /2. In this case, we must replace the cosine func- 
tion of Eq. (jBlj) with a sine function, owing to the con- 
jugation of quantized currents and voltages. The max- 
imum vacuum voltage of, e.g., resonator A is given by 
v A0 = y / hoj A /2C 1 . li ~ 0.5871 U-V for the realistic param- 
eters of Tabic HI At T^sim/2, we then obtain the max- 
imum voltages in the simulated regions, which are ap- 
proximately =f0.065v ao . Towards the center of the simu- 
lated regions the voltage reduces to zero because its cor- 
responding first mode is characterized by a node. Also, 
higher modes do not contribute for the same detuning ar- 
guments outlined above. Therefore, we can safely neglect 
all capacitive couplings in our simulations. 

Second, we notice that £ sim = 500 |J.m is chosen to be 
large enough compared to the lateral dimension Z ql = 
200 urn of the flux qubit loop [see Figs. EUa)-(c)]. This 
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avoids errors due to fringing effects when simulating the 
coupling coefficients between qubit and resonators. For 
consistency, we have also performed several simulations 
assuming larger values of ^ sim , up to 1-1.5 mm (data not 
shown) j££ We have not found any appreciable deviation 
in the resulting inductances. 

Third, we stress that special care has to be taken when 
using FASTHENRY to simulate the second-order induc- 
tances of our three-circuit network. In order to compute 
the inductance matrix, test currents must be applied to 
the involved structures at specific probing ports. How- 
ever, when applying test currents to all three circuits si- 
multaneously, only first-order inductances are calculated. 
This is due to the boundary conditions that must be ful- 
filled at the probing ports. This fact has important im- 



plications for the calculation of the mutual inductance 
to, which is the sum of first- and second-order mutual 
inductances between resonators A and B. In this case, 
it is crucial to apply test currents only to the two res- 
onators, but not to the qubit circuit. On the contrary, 
the pure first-order mutual inductance m between A and 
B can be simulated in two equivalent ways: Either the 
qubit circuit is completely removed from the network or 
test currents are applied to all three structures. We do 
not notice any difference between these two approaches. 
The above arguments also apply to the calculation of the 
renormalized self-inductances L* a and L* h of the two res- 
onators and their pure counterparts L*. a and L* h . 
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